LOWER BOUND ON THE DENSITY OF STATES FOR 
PERIODIC SCHRODINGER OPERATORS 



SERGEY MOROZOV, LEONID PARNOVSKI, AND IRINA PCHELINTSEVA 

Abstract. We consider Schrodinger operators — A + V in K d 
(d ^ 2) with smooth periodic potentials V and prove a uniform 
lower bound on the density of states for large values of the spec- 
tral parameter. 



1. Introduction 

Let H = — A + V be a Schrodinger operator in L2(M. d ) with a smooth 
periodic potential V. We will assume throughout that d ^ 2. The 
integrated density of states (IDS) for H is defined as 

N(X) := lim L~ d N(X;H^), XeR. (1.1) 

L^oo 

Here Hp is the restriction of H to the cube [0, L] d with the Dirichlet 
boundary conditions, and iV(A; ■) is the counting function of the dis- 
crete spectrum below A. For H := — A the IDS can be easily computed 
explicitly (e.g. using the representation ( 12.61) below): 

iV (A) = (2n)- d d- W d/2 - (1-2) 

Here Ud = 2ix d l 2 /T(d/2) is the surface area of the unit sphere S^" 1 in 
R d . 

The asymptotic behaviour of the function f 1 1.1 1) for large values of the 
spectral parameter was recently studied in a number of publications, 
see PP, [3], and references therein. 

Our article concerns the high-energy behaviour of the Radon-Niko- 
dym derivative of IDS 

g := dN/dX, 

which is called the density of states (DOS) (see jl]). Our main result 
is that for big values of A 

g(X) > 9o(X) (1- o(l)), (1.3) 

where 

g (X) = dN (X)/dX = (27rr<W- 2)/ 72. 
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We remark that (11.31) should be understood in the sense of measures; 
in particular, we do not claim that g(X) is everywhere different iable. 

It has been proved in [2j that the spectrum of H contains a semi-axis 
[Ao, +oo); this statement is known as the Bethe-Sommerfeld conjecture 
(see the references in [2] for the history of this problem). This result 
has an obvious reformulation in terms of IDS: each point A ^ Ao is a 
point of growth of N. It was also proved in [2] that for each neN and 
e = A~ n we have 

N(X + e) - N(X) < e\ id - 2)/2 . (1.4) 

Later, when the second author discussed the results and methods of 
[2] with Yu. Karpeshina, she suggested that using the technique from 
that paper, one should be able to prove the opposite bound 

N(X + e) - N(X) > eX {d - 2)/2 (1.5) 

when A is sufficiently large, not just with e = X~ n (when the proof 
is relatively straightforward given [2]), but also uniformly over all e G 
(0, 1]. In this paper we prove that for big A 

N(X + e)- N(X) * ^y^ (d - 2)/2 i 1 ~ • (L6) 

Note that (jl.6p implies the claimed bound (11.31) . 

The proof of (jl.6p is heavily based on the technique of [2] and uses 
various statements proved therein. In order to minimise the size of our 
paper, we will try to quote as many results as we can from [2] , possibly 
with some minor modifications when necessary. 

Acknowledgement. As we have already mentioned, this paper 
is a result of observations and suggestions made by Yu. Karpeshina; 
we are very grateful to her for sharing them with us and allowing 
us to use them. The authors were supported by the EPSRC grant 
EP/F029721/1. 

2. Preliminaries 

We study the Schrodinger operator 

H = -A + V(x), xGR d (2.1) 

with the potential V being infinitely smooth and periodic with the 
lattice of periods A. We denote the lattice dual to A by A', fundamental 
cells of these lattices are denoted by fl and fl\ respectively. We choose 
to be the first Brillouin zone and introduce 

Q :— sup |£ G (2.2) 
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Let 

D := -iV, D(k):=D + k. (2.3) 
The Floquet-Bloch decomposition allows to represent our operator 
direct integral (see e.g. [1]): 

H= I ©#(k)dk, (2.4) 

where 

tf(k) = D(k) 2 + V(x) (2.5) 
is the family of 'fibre' operators acting in L 2 (fi). The domain of each 
H(k) is the set of periodic functions from H 2 (Q). The spectrum of H 
is the union over k G ft* of the spectra of the operators (12.51) . 

We denote by | • | Q the surface area Lebesgue measure on the unit 
sphere S^ 1 in R d and put u d := \S d -% = 2n d ^ 2 /T(d/2). Finally, 

N(X) := (2n)- d [ #{j : A,(k) < A} dk (2.6) 

is the integrated density of states of the operator (12. ip . It is known (see 
e.g. [4J) that the definitions (II. ip and (I2.6P are equivalent. 
The main result of the paper is 

Theorem 2.1. For sufficiently big A and any e > the integrated 
density of states of H satisfies (II. 6p . 

By B(R) we denote the ball of radius R centered at the origin. Given 
two positive functions / and g, we say that / g, or g <C f, or 
g = 0(f) if the ratio gj f is bounded. We say / x g if / g and /Cg. 
Whenever we use O, o, -C, or x notation, the constants involved can 
depend on d and norms of the potential in various Sobolev spaces H s ; 
the same is also the case when we use the expression 'sufficiently large'. 
By A = p 2 we denote a point on the spectral axis. We also denote by v 
the Loo-norm of the potential V, and put J := [A — 20v, A + 20v}. Let 

A := {£ G M d , 1 1£| 2 - A| < 4(h;}. (2.7) 

Notice that the definition of A obviously implies that if £ G A, then 

Any vector £ G M d can be uniquely decomposed as £ = n + k with 
n G A^ and k G fiL We call n = [£] the 'integer part' of £ and k = {£} 
the 'fractional part' of 

By vol(-) we denote the Lebesgue measure in Mr. The identity matrix 
is denoted by I. For any h G L^iS^) we introduce its Fourier coefficients 

h n := (vo\ny 1/2 / /i(x)exp ( -i(n,x)) dx, n G A f . (2.8) 
Jn 
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For £ G M. d \ {0} we define r = r(£) := |£| and g := We put 

R = R{p) : = p V(36^(d+2)) (2 g) 

(so that the condition stated after equation (5.15) in [2 J is satisfied). 
For j eN let 

e; := At fl B(jR) \ {0}. 
Let M := 5d 2 + 7d. We introduce the set 

2 := G .A| | | > P 1/2 , for all r/ G 0' 6M }. (2.10) 

In other words, !B consists of all points £ G .A the projections of which 
to the directions of all vectors rj G 0' 6M have lengths larger than p 1 / 2 . 
We also denote D := A \ 'B. 

In the rest of the section we quote some results from [2] which we 
will use in this paper. Our approach is slightly different from that 
of [2]. In particular, we consider arbitrary lattice of periods A, not 
equal to (2vrZ) d . We also use a different form of the Floquet-Bloch 
decomposition (so that the operators on fibers (12. 5p are defined on the 
same domain). This leads to several straightforward changes in the 
formulation of the results from [2]. These changes are: 

(1) The lattices (2irZ) d and Z d are replaced by A and At, respec- 
tively. The 'integer' and 'fractional' parts are now defined with 
respect to At (see above); 

(2) The matrices F and G are replaced by the unit matrix I through- 
out; 

(3) The Fourier transform is now defined by (12.81) . and the expo- 
nentials e m introduced at the beginning of Section 5 in [2] are 
redefined as 

e m (x) := (volfi)- 1/2 e i(m ' x) , m G A f ; 

(4) The operators H(\t) are now given by (12.51) on the common 
domain D. 

The main result we will need follows from Corollary 7.15 of [2]: 

Proposition 2.2. There exist mappings f,g:A—>M. which satisfy 
the following properties: 

(i) /(£) is an eigenvalue of H(k) with {£} = k; — |£| 2 | ^ 2v. 
f is an injection ( if we count all eigenvalues with multiplicities ) and 
all eigenvalues of H(\s) inside J are in the image of f . 

(ii) If £eA, then \ f{£) - g{£)\ ^ p~ d ~ 3 . 
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(Hi) For any £ G S 

a(i) = III 2 

+E E E c ». ,<«.^>-«...<e,^. (!! - 11) 

J'=l Tj 1 ,...,T 7j ee' A/ 2<Cm+-+n 3 <2M 

('raj Let J = [a, b] C A be a straight interval of length L : = |b — a| <C 
p^ 1 . Then there exists an integer vector n such that gr(b+n)— gr(a)| <C 
Lp + p~ d ~ 3 . Moreover, suppose m ^ is an integer vector such that the 
interval I + m is also entirely inside A. Then there exist two different 
integer vectors rix and n 2 such that \g(h + ni) — g(sL)\ <C Lp + p~ d ~ 3 
and \g(h + n 2 ) — g(a + m) I <C Lp + p~ 



-d-3 



Remark 2.3. Formula (12.111) implies that 

dg/dr(£) x p, for any £ e B. (2.12) 

For each positive 5 ^ t> we denote by A(S), B(<5), and 2) (5) the 
intersections of g _1 ([p 2 — 5, p 2 + 5}) with .A, 25, and D, respectively. 
It is proved in Lemma 8.1 of [2 J that 

vol (D{6)) < p d - 7/3 S. (2.13) 

The following statement (Corollary 8.5 of [2]) gives a sufficient con- 
dition for the continuity of /: 

Lemma 2.4. There is a constant C\ with the following properties. Let 

I : = {£(t) : t G [t mi n,t max }} C 

fre a straight interval of length L < p~ 1 5. Suppose that there is a 
point to G \t m in,t max ) with the property that for each non-zero n G A" 1 " 
<?(£(to) + n) zs either outside the interval 

'g(t(t )) - c lP - d - 3 - c iP l, g (t(t Q )) + c lP - d - 3 + C lP L 

or not defined. Then /(£(£)) a continuous function oft. 

By inspection of the proof of Lemma 8.3 of [2] we obtain 

Lemma 2.5. For Zetrye enough p and 5 < p _1 i/ie following estimates 
hold uniformly over a G A" 1 " \ {0}: if d ^ 3, 

vol n + a)) < (5 2 p d_3 + <*P~ d ); (2.14) 
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*fd = 2, 

f<<5 3 / 2 , |a|<2p-l, 
vol(3(<5)n (3(5)+ a) W <5 3 / 2 + 5p- 2 , ||a|-2p|<l, (2.15) 

[=0, |a|^2p+l. 

3. Prevalence of regular directions 
Lemma 3.1. For p big enough and 

< 5 «C /T^ 3 

i/iere exists a set 3 = ?(p) on the unit sphere S^ 1 in M. d with 

|5F| ^ d (l-o(l)) (3.1) 

such that /(£) is a simple eigenvalue o/ //"({£}) continuously depending 
onr := |£| for every £ = (r, £') G f' 1 ([p 2 -S, p 2 + 5)) with £' : = G 
J. 

Proof. It is enough to consider 5 := p~ d ~ 3 . For each £' G S^ -1 let 

%W:=K,r>0}nS(5). (3.2) 

Let ^ := {£' G S^/^) ^ 0, p)nD((5) = 0}. 

For any t] G 0' 6M the area of the set of points £' G S d ~ l satisfying 



|<r*W)| < 



1/2 



is evidently 0(p~ 1 ^ 2 ) if r ^ p/2 (the latter is true for all r£' G A). 
Since the number of elements in Q' eM is 0(R d ), by (12. 9p and (12.101) we 
have 

l^- 1 \^x| = o(l). (3.3) 

By definition 3(5) = 3 n ^([p 2 - 5, p 2 + 5]), hence ( 12+2|) implies 
that for big p the length lg(S) of I^'(5) satisfies 

l $/ (5)^S P ~\ C'eJi. (3.4) 

Let 

3" := {£' G 3^ | / is continuous on Ig(6)}, 

and 

Lemma 12.41 tells us that for each point £ G £(5) there is a non-zero 
vector n G A" 1 " such that 

\9^ + n)-g($)\^C 1 (p- d - 3 + pl $ i5)) «(p- d - 3 + <5). (3.5) 

Since \g(£) — p 2 \ ^5, this implies 

\g{£ + n) - p 2 | < C 2 (p- d " 3 + 5) =: 5r « p^" 3 = 5, 
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and thus £ + n G A(5i); notice that C2 > 1 and so Si > 5. Therefore, 
each point £ E 8(5) also belongs to the set (A(Si) — n) for a non-zero 
n G A" 1 "; obviously, |n| <C p. In other words, 

8,(5) C |J (A(5 1 ) - n) = |J (30*0 - n) U |J (©(50 - n). 

neAtnB(Cp),n^O n^O n^O 

(3-6) 

To proceed further, we need more notation. Denote Dq(Si) to be the 
set of all points v from D(5i) for which there is no non-zero n G At 
satisfying v — n G 3(5); Di(5i) to be the set of all points v from D(5\) 
for which there is a unique non-zero n G A^ satisfying u— n G S(<5); and 
1)2(^1) to be the rest of the points from D(5i) (i.e. D 2 (5i) consists of 
all points v from D(5i) for which there exist at least two different non- 
zero vectors n!,n 2 G A^ satisfying v> — rij G 23(5)). Then Lemma 8.7 
of [2] implies that we can rewrite (13.61) as 

8(5) C |J (B(5i) ~ n) U |J (Di(5i) - n). (3.7) 

n^O n^O 

This, obviously, implies 
8(8) C |J ((^)-n)nS(5))uU ((dCtfO -n)n 3(5)), (3.8) 

n^O n^O 

since 8(5) CB(5). 

The definition of the set T>i(8i) and (I2.13P imply that 




«C vol (D(<Ji)) < 5ip d - 7/3 < (5p d ~ 7/3 . 

For <i ^ 3 Lemma I2.5[ inequality 5 < 5\ , and the fact that the union 
in (13.81) consists of no more than Cp d terms imply 

vol ( |J ((3(50 - n) n ^(Sjj) < /(^V- 3 + 5 ip - d ) < 5(p d ~ 6 + 1). 

(3.10) 
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For d = 2 we obtain by Lemma 12.51 

vol( |J (s(5)n(S(5!) + n) 



neAt\{0} 



^ vol(s(5)n (3(5i) + n 



neAt\ {0 } 

M<2p-i (3.11) 



neAt\{0} 
||n|-2p|<l 

<<£V + p(^ + *ip- a )<*p- 1/a , 

where we have used that 

#|n G A 1 " ||n| - 2p| < l| < p. 

Applying M . fl3~T0l) . and (EH]) to (ESD we obtain for all d ^ 2 

vol £(<5) <5/- 7/3 . (3.12) 

By definition, 

£(*)= IJ W 

Hence by ((321) 

|ff' 1 \ff'| <rV- ,, vol£(<J). (3.13) 
Combining (13.121) and (13.131) we conclude that for big p 

|2i\31o = o(l). (3.14) 

We have 

IS*' 1 \ ?\ = IS"*" 1 \ 5 X | Q + |?i \ 5| Q . (3.15) 

Substituting (ED and (BHD into (EED we obtain (ETIli . 

Now we notice that for every ^' 6 J the interval Ig(5) has the 
following property: for each point £ G (5) and each non-zero vector 
n G At such that £ + n G .A we have |g(£ + 

n ) > 2p-^ 3 . This 

implies /(£ + n) — /(£) ^ 0. Therefore, /(£) is a simple eigenvalue of 
for each £ G ^'(5). The lemma is proved. □ 

4. Some properties of operators on the fibers 
For m G M let 

y(m) ._ f^J2 \n\ 2m \V n \^ ll J ' 

neAt 
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Since V is smooth, is finite for any m ^ 0. Recall that Q is 

defined by Q. 

Lemma 4.1. Fix m G N and x G (0,1). For k G Zei ip be a 

normalized eigenf unction of H(k): 

H(k)i/; = (4.1) 

with the eigenvalue 

C > max{36g 2 x- 2 , (l + mx)^^*}. (4.2) 
Then there exists M m = M m (d,A,V) G R+ swc/i that for all n G 

|n| ^ (l + mx)V7 (4-3) 
t/ie Fourier coefficients oft]) satisfy 

|V»n| < M m x- m |n|- (3m+1)/2 . (4.4) 

Proof. We proceed by induction. Suppose that either m = 1, or m > 1 
and the statement is proved for m — 1. Substituting the Fourier series 

V>(x) = (volfi)~ 1/2 ^ ip n exp (i(n,x)), x£!l 

neAt 

into (14. 1 p and equating the coefficients at exp (i(n, x)J on both sides, 
we obtain by (12.51) : 

|n + k|>n + X) y n-l^l = C^n. (4.5) 

leAt 

Since |k| < Q, by (0~2D and (USD we have 

2|n||k| «C x|n| 2 /6 + Qx^Q 2 ^ x|n| 2 /3. (4.6) 
For x G (0, 1), it follows from flO]) that 

|n| 2 -C ^ (l-(l + x)" 2 )|n| 2 = x(2 + x)(l + x)~ 2 |n| 2 ^ x|n| 2 /2. (4.7) 
Combining (I4.6[) and (14.71) we obtain 

|n + k| 2 -( > |n| 2 - 2|n||k| - £ > x|n| 2 /6, 
and thus by (14. 5 p 

|V„| <6x- 1 |n|- 2 ^|K-i^i|. (4.8) 
leAt 

If m — 1 we estimate the sum on the r. h. s. by using Cauchy- 

Schwarz inequality (since ip is normalized) and obtain (14.41) with Mi : = 
6 y(o). 
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If m > 1, we estimate 

£ IK-iViK sup |^ m | £ (4.9) 

leAt:H-n|<|n|Vi m: |m|^|n|-|n|l/«« leA t : |l|<|n|V«i 

By (14.31) . (I4.2p . and monotonicity of the function q(t) = t — t l l d for 
t > 1 we have 

|n| - (nl 1 ^ ^ (l+mx)^ - ((l + mx)^) 1 ^ (l + (m - l)x) yft. 
We can thus apply the induction hypothesis obtaining 



SUp |V> m | ^ X^M^l - lnjCl-^/") | n 



m:|m|^|n|-|n|l/d 



(1-d) /d^ l-3m/2 | ^ | l_3 m /2 

(4.10) 



Since x £ (0, 1), by QOJ and we have 
|n| ^ (1 + mx 

thus 

(1 - |n|( 1 - d )/ d ) 1 - 3m/2 < 2 3 " 1 / 2 - 1 . (4.11) 

Let 

W :=supr~ d #{\e A + | |1| < r}. 
Clearly, < oo. Then by Cauchy-Schwarz inequality 

M^W^V^^ 2 . (4.12) 

leAt: |lK|n|Vd 

Substituting (OOj) . and (14TT2]) into (jMD we get 

l^-'^l < 2 3m/2 - 1 x 1 - m ^ 1 / 2 \/WM m _ 1 |n| 3 ( 1 - m )/ 2 . 

leAtiH-nKlnlVd 

(4.13) 

On the other hand, since ||^|| = 1, by Cauchy-Schwarz we have 

IK-iViI < In^-)/ 2 £ lll^^lViH^-il 

leAt: |I-n|>[n|Vd leAt: |l|>|n|!/ d 

^ | n |3(l-m)/2^ ^ 1 1 j 3(m— l)d! | | 2 ^ < ^(3(m-l)d/2) | n |3(l-m)/2_ 

leAt: |l|>|n| 1 / d 

(4.14) 

Inserting (14.131) and (14. 141) into ( l4~8l) we obtain with 
M m := 6(2 3m / 2 - 1 iy 1 / 2 y^M m _ 1 + yWm-Dd/a)). 

□ 



|V k C(k)| < 2||D(k)V(k)|| = 2( |n + k| 2 |^ n (k)| 2 ) V2 . (4.18) 
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Lemma 4.2. For any rj e (0, 1) there exists Co > such that if ((h) 
Co is a simple eigenvalue of H(k) for some kefit then 

|V k C| < 2(1 + 77)^/?- (4.15) 
Proof. Let i/>(k) be the eigenfunction corresponding to C(k) with 

|^(k)|| = l. (4.16) 

Then 

V k C(k) = V k (^(k),F(k)^(k)) = (>(k), (V k ff(k))V>(k)). (4.17) 
By (E3D and Q, 

V k -ff(k) = 2D(k). 
Substituting this into (14.171) we obtain: 

|n + kr^ n (k)|-J 1/ ' 

ieAt 
Let 

m := [(d + l)/3] + 1 (4.19) 

and 

x:=r//(2m + l). (4.20) 

We put 

Co := max {36Q 2 x- 2 , (1 + mxP^i)) ( 4 .21) 

and assume that C := C(k) > Co- Since by (E2j) |k| ^ Q, by (jOBjl . 
ffl~211 . and (14^201 we have 

J2 |n + k| 2 |^ n (k)| 2 < (l + (m + l/6)x) 2 C< (l + r//2) 2 C- 

|n|<(l+mx)V? 

(4.22) 

For |n| ^ (1 + mx)\/C we apply Lemma [4.11 obtaining 

J2 |n + k| 2 |^ n (k)| 2 < M 2 m x~ 2m J2 \n + k\ 2 \n\- 3m -\ 

|n|>(l+m*)V£ |n|>(l+mW)-v/f 

(4.23) 

By (14.191) the r. h. s. of (14.231) is finite and is 0(C -1 ^ 2 )- Thus, increasing 
Co if necessary, by (Qgj) . (14^221) . and (OBI we obtain (1435]) . □ 
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5. Proof of Theorem 12.11 

It is enough to prove 

Theorem 5.1. For any a G (0, 1) there exists po > big enough such 
that for all p ^ p 

N(p 2 + 8)- N(p 2 - 8) > (1 - a)(2it)- d iu d 6p d ~ 2 (5.1) 

for any 

< 8 ^ p- d ^. (5.2) 

Indeed, the original statement of Theorem 12.11 can be obtained by 
partitioning of the interval [A, A + e] into subintervals with lengths not 
exceeding 2\(~ d ~ 3 ^ 2 and adding up estimates ( 15.11) on this subintervals 
(with p 2 being respective middle points). 

Proof. We first express the growth of IDS in terms of the function / of 
Proposition 12.21 

N{p 2 + 8)- N{p 2 -8) = (27r)- d vol (rV - S,p 2 + 8)). (5.3) 
We can write 

r poo 

vol (r 1 [ p "-8,p 2 + 8))= / X (r, £ V" 1 dr d£', (5.4) 

Js^ 1 Jo 

where x is the indicator function of / _1 ([p 2 — 8,p 2 + 8)). To obtain a 
lower bound we can restrict the integration in ( 15.41) to £' G 5" defined 
in Lemma I3.ll Then for any rj G (0, 1) there exists po > such that 
for any p ^ p we have 

|5| > (l-v)^d, (5.5) 

and for any £' G 5F the support of x(">£') contains an interval [ri,r 2 ] 
with 

(1 — 7^)p ^ r! < rx + (1 — Ti)p~ l 8 ^ r 2 . (5.6) 

Indeed, the first inequality in (15.61) follows from Proposition 12.2( h) . (iii) . 
The last inequality in (15.61) follows from Lemmata 13.11 and 14.21 
Thus for all p ^ p by (15.51) and (15.61) we obtain 

Xir^y-Ur&l-' > [ (l- V ) d p d - 2 8d£ 

J 3 (5.7) 

>(l-r)) d+1 u; d p d - 2 8, 

Combining ( 15. 3L ( 15. 4p . and ( 15.71) . and choosing r\ small enough we 
arrive at (15.11) . The theorem is proved. □ 
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